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The magnetic phase diagram and the thermodynamic properties of a square lattice containing 
antiferromagnetically coupled spin-1/2 triangles are studied. A Heisenberg Hamiltonian with three 
strong intra-triangle exchange coupling constants and one weak inter-triangle exchange coupling 
constant is proposed to model the interacting Cu 2+ ions of the CU3M0O4 planes in La4Cu3MoOi2. 
Depending on the ratio of the intra-triangle coupling constants, various long-ranged magnetic phases 
are shown to compete at low temperatures. A comparison of numerical calculations of thermody- 
namic properties in the model Hamiltonian with recent experiments on La4Cu3MoOi2 suggests that 
the spins in this material order antiferromagnetically along the x-direction and ferromagnetically 
\f~) • along the y-direction of the CU3M0O4 planes. The temperature dependence of the uniform magnetic 

susceptibility, the heat capacity, and the entropy are calculated and are shown to be in qualitative 
agreement with the experiments on this compound. 

It is known that the zero-point fluctuations in two-dimensional (2D) quantum antiferromagnets are responsible for 
a strong renormalization of their order parameter., causing a substantial correction of the sublattice magnetization 
in long-ranged antiferromagnets such as La2Cu04U, or completely suppressing the magnetic moment in short-range- 
ordered antiferromagiietic (AF) compounds such as the structurally frustrated Kagome lattice, containing coupled 
spin-1/2 trianglesa □ In contrast to systems with long-range AF order, this spin liquid is best described within a 
£h ' resonant valence bond picture that accounts for the degeneracy of the ground state, the singlet-triplet spin gap, and 
the large number of low-lying singlet states in the thermodynamic limitBij. n 

i^i - Recently, thermodynamic measurements on the triangular cluster compound La4Cu3MoOi2 were reportedu. This 

system is a layered material, containing planar CU3M0O4 square lattices of nearly isosceles Cu 2+ spin-1/2 triangles, 

*~H . as shown in Fig. 1. In contrast to the Kagome antiferromagnet, La4Cu3MoOi2 orders antiferromagnetically at a Neel 
temperature of 2.6K despite the structural frustration within the triangles. This indicates that - unlike the Kagome 
lattice - the bipartite nature of the square lattice of triangles in this compound is essential for the long-range ordering. 
Associated with the antiferromagnetic ordering transition, an entropy change of 1.56 J/(K mol Cu) was observed. 
Furthermore, a plateau in the magnetization above a critical field of approximately 20 T was reported. 
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FIG. 1. Square lattice of Cu 2+ triangles in the CU3M0O4 planes of La4Cu3MoOi2. The intra-triangle couplings Ji, J2, and 
J3 are strong compared to the inter-triangle couplings J'. 



From the different Cu-0 bond lengths within and between the triangles|3 it can be concluded that the intra 
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triangle couplings J\, J2, and J3 are much stronger than the inter-triangle coupling J'. The magnetic moments of 
the neighboring Cu 2+ ions within the CU3M0O4 planes interact via antiferromagnetic superexchange, mediated by 
the filled O 2- p-orbitals. This leads us to propose a planar model Hamiltonian of weakly coupled spin-1/2 trimers 
for La4Cu3MoOi2, given by 



H — 2_^ [JlSr,l • S r _2 + J2&r,2 • S r ^3 + J3S r: 3 • S r> i + J (S r ,3 • S r 



-x,l 



r,2 ' S r +yj 



r,2 • S r+ y^ 3 ) 



(1) 



where S rj - represents a spin-1/2 degree of freedom at site j of the triangle centered at position r. The underlying 
2D square lattice of triangles is spanned by the unit vectors x and y. In order to account for deviations from perfect 
isosceles triangles, 3 different intra-triangle couplings-,are considered. From experience, it is rather difficult to derive 
the precise exchange constants from first principles.B An alternative approach is to determine them by comparing 
results of model calculations on H to experimental data. However, this approach is typically complicated by the 
fact that various parameter sets may fit the experimental results equally well.Q In the following, we will therefore 
first examine the weak-coupling phase diagram of H. Subsequently, the recent experiments on La4Cu3MoOi2 will be 
analyzed, based on numerical calculations for the proposed model Hamiltonian in the relevant parameter regime. 

The weak-coupling (J' <C Ji, J2, J3) zero-temperature phase diagram of H is determined by the dominant spin-spin 
correlation function of the total trimer spins S r = S r i + S r 2 + S r 3, defined by 5 q = X) r ex P (*1 ' r )(So ■ S r ). We 
have calculated this quantity by exact numerical diagonalization of finite clusters with periodic boundary conditions 
for small inter-triangle coupling strengths J' .a It is found that the wave vector of the dominant spin-spin correlation 
function strongly depends on the relative strengths of the intra-triangle couplings J\, J2, and J3. The resulting zero- 
temperature phase diagrams for square lattices of weakly coupled triangles are shown in Fig. 2 with (a) J' = O.OIJ3, 
and (b) J' — O.IJ3, where the various phases are labeled by their magnetic ordering wave vectors Q = (q x ,q y ). 
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FIG. 2. Zero-temperature magnetic phase diagram of the AF spin-1/2 trimer square lattice with a weak inter-triangle 
coupling (a) J' = O.OU3 and (b) J' = O.U3. 



Since the model Hamiltonian H is invariant under J\ <-» J2, these phase diagrams have reflection symmetry about 
the diagonal axis J\ = Ji- AH exchange coupling constants have been scaled by the intra-triangle coupling constant 
along the x-direction, J3. In Fig 2(a), the phase diagram is shown for the ultra-weak-coupling limit with a small inter- 
triangle coupling J' = O.OU3. At the isotropic point J\ = J2 = J3, the system orders antiferromagnetically in both 
planar directions, resulting in an ordering wave vector Q = (tt, 7f). For small J3 (J3 < J\, J2), the system decouples 
into spin-1/2 AF chains along the x-direction with weak residual ferromagnetic (FM) couplings along the y-direction, 
leading to Q — (tt, 0). In contrast, in the regimes J\ 3> J2, -h and J2 3> Ji, J3, the model describes effective spin-1/2 
AF chains along the y-axis with weak residual FM coupling along the x-direction, and thus Q = (0, tt). The symmetry 
line J\ — J2 with Q = (tt, 0) is found to separate these two Q = (0, tt) phases.Q 

As the inter-triangle coupling J' is increased (Fig. 2(b)), a finite region with Q = (jr, tt) opens up, centered around 
the isotropic point J\ = J2 = Jz- Furthermore, the Q = (tt, 0) regime is found to extend away from the diagonal 
J\ = J2 line. It has been proposed that the AF ordered low-temperature phase of La4Cu3MoOi2 may fall into the 
Q = (tt, tt) regimejj However, as we will discuss in the following, a numerical analysis of the thermodynamic properties 
of H - in particular the magnetic entropy - suggests that this compound should be placed in the regime of weak J3 
with Q = (tt, 0) rather than in the Q = (tt, tt) regime close to the isotropic point. 

When analyzing the thermodynamic response functions for the square lattice of weakly coupled AF spin-1/2 trian- 
gles, three temperature regimes can be identified. At high temperatures T > J±, J2, J3, the system is a paramagnetic 



ensemble of thermally decoupled spins. Upon lowering the temperature into the region J 1 < T < Ji, J2, J3, the intra- 
triangle couplings become relevant, and spin doublets are formed on each triangle. The properties of the system in 
this temperature interval are thus well characterized by a thermodynamic ensemble of decoupled trimers with effective 
spin-1/2 degrees of freedom. Below the ordering transition temperature TV « J' the small inter-triangle couplings 
correlate these spin-1/2 trimers, leading to long-range-ordered magnetic phases with characteristic Q vectors that 
depend on the ratio of the intra-triangle coupling constants, as discussed above. We have determined the temperature 
dependence of the uniform magnetic susceptibility, the heat capacity, and the entropy from full numerical diagonal- 
izations of H on finite clusters with periodic boundary conditions. Results representative for the the Q = (tt, 0) phase 
with (Ji, J 2 , J3) = (1.5,1.5,1.0), Q = (0, 7r) phase with (Ji, J2, ^3) = (0.5,1.0,1.0), and the Q = (7r, 7r) phase with 
(Ji, J2, J3) = (1.0, 1.0, 1.0) are shown in Fig. 3. 



0.5 



0.4 



0.3 



h- 



0.2 



0.1 



0.0 



0.10 




,^ 


0.05 


\ 


- 


0.00 
C 






) 


0.2 0.4 


tt S 


't^^' 




(a) 



0.0 




CO 



1.0 2.0 3.0 4.0 

T/J, 



(b) 







0.4 


/ 
/ 
/ 
/ 

/ 




0.2 


¥'/^~ 


- 




1 / 

1/ 


(c) 



0.0 1.0 2.0 3.0 4.0 
T/J' 



FIG. 3. Thermodynamic response in the square lattice of AF spin-1/2 triangles, (a) Temperature-dependent uniform sus- 
ceptibility in the Q = (7r,0) regime (solid line), the Q = (Q,7r) regime (dot-dashed line), and the Q = (it, it) regime (dashed 
line). The inset shows the intermediate-temperature plateaus of yT at j^. (b) Temperature dependence of the specific heat, 
and (c) the magnetic entropy per site. 



The three distinct temperature regimes are clearly observed in the uniform susceptibilities, shown in Fig. 3(a). 
At high temperatures the spins are effectively uncorrelated with a Curie constant lim-r^oo [%T] = C = i- At 
intermediate temperatures there is a plateau in %T with a reduced effective Curie constant C" = yj, associated with 
the paramagnetic response of uncorrelated spin-1/2 trimers, containing three spins each. Finally, at low temperatures 
of order J' these trimers condense into magnetically ordered phases, leading to an additional reduction in \F . These 
trends in the temperature dependence of the uniform susceptibilities are similar in all three parameter regimes. 
Furthermore, the shape of x(T) is found to be in good qualitative agreement with measurements on La4Cu3MoOi2 
(see Fig. 2 of Ref.H), which show an onset of magnetic ordering at 2.6 K, two orders of magnitude smaller than the 
observed threshold temperature for trimer formation at approximately 250 K. 

In Figs. 3 (b) and (c) the specific heat and the magnetic entropy per site are shown. The onset of the intermediate 
regime in which intra-triangle correlations become relevant is marked by the maxima in the specific heat divided by 
the temperature, C/T, consistent with the onsets of the plateaus in xT, seen in Fig. 3(a). For the Q = (tt, 0) and 
the Q = (0, tt) regimes, the entropies per site exhibit plateaus at S = | In 2, rising up to In 2 at high temperatures 
(not shown in the figure). Hence, only 1/3 of the magnetic entropy is released in the intermediate temperature 
regime, corresponding to an ensemble of N/3 trimer doublets with a decreased paramagnetic moment of spin-1/2 per 
trimer. In contrast, the entropy for the isotropic coupling case with Q = (tt, tt) shows a larger plateau value of , In 4, 
associated with the degeneracy of two Kramers doublets in this regime. In fact, the same plateau value is found over 
the whole Q = (it, tt) regime centered around the isotropic coupling point. At very high temperatures the entropy 
saturates at the. fully paramagnetic value of In 2 in all cases. A comparison with the experiments on La4Cu3MoOi2 
(Fig. 3 in RefJj) clearly favors the smaller intermediate-temperature plateau value of 5 In 2 over ^m4, suggesting 
that the low-temperature magnetic order in this material is characterized by Q = (tt, 0) or Q = (0, tt) rather than by 

Q = (7T,7r). 

Examining the intra-triangle bond lengths and bond angles in the CU3M0O4 planes of La4Cu3MoOi2 (table 1 of 



Ref.Ef), it appears likely that J3 is smaller than both J\ and J2BE2I This suggests that the magnetically ordered 
low-temperature phase in this compound is part of the Q = (n, 0) regime rather than the Q = (0, 7r) regime. Neutron 
scattering studies should help to verify this conjecture. 

In summary, we have studied the magnetic phase diagram and the thermodynamic response of a square lattice 
containing weakly coupled AF spin-1/2 triangles. A Heisenberg Hamiltonian with three strong intra-triangle exchange 
coupling constants and a weak inter-triangle exchange coupling constant is proposed to model the interacting Cu 2+ 
ions of the CU3M0O4 planes in La4Cu3MoOi2. Depending on the ratio of the intra-triangle coupling constants, 
various long-ranged magnetic phases are possible at low temperatures. A comparison of numerical calculations of 
the thermodynamic properties in the proposed model Hamiltonian with recent experiments on La4Cu3MoOi2 leads 
us to suggest that the magnetically ordered phase in this compound belongs to the Q = (n, 0) regime rather than 
the Q = (0, 7r) or the Q = (ir, it) regimes of the model. The ordering wave vector could be determined in a neutron 
scattering study, giving conclusive information regarding the nature of the magnetically ordered phase in this material. 
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